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During this period, most of the issues pertaining to the Hyperbolic 
Navier-Stokes equations have been sorted out and an Abstract (copies en- 
closed) was submitted to the CFD meeting in Honolulu. 

Work is continuing on incorporating the viscous terms into the re- 
acting flow code. When this is completed and checked, we intend to in- 
corporate a chemistry model based on the complete react ^ on# More- 

over, all transport properties will be based on the 12-6 potential. 
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ABSTRACT 

Using ideas fro* kinetic theory, the Navier-Stokes equations are modified 
in such a way that they can be cast as a set of first order hyperbolic 
equations. This is achieved by incorporating time dependent terms into the 
definition of the stress tensor and the heat flux vectors. The boundary 
conditions are then determined from the theory of characteristics. Because 
the resulting equations reduce to the traditional Navier-Stokes equations when 
the steady state is reached, the present approach provides a straightforward 
scheme for the determination of inflow and outflow boundary conditions. The 
method is validated by comparing its predictions with known exact solutions of 
the steady Navier-Stokes equations. 
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Introduction and Approach 


Solutions of the Navier-Stokes equations have not matured to the point 
where one can routinely simulate transition, separation and turbulence. One 
of the major obstacles that prevents us from achieving these goals is our 
inability to formulate accurate inflow and outflow boundary conditions . 1 
Attempts at formulating such conditions centered around the assumption that, 
in the far field, the governing equations reduce to a simple form. Such forms 
may include a wave equation, a Poisson's equation or, some other simple 
equation. Asymptotic solutions of these equations are then constructed and 
used to provide boundary conditions for the Navier-Stokes equations at the 
computational boundaries . 2-0 

Because viscous effects are important in the vicinity of bodies and in 
their wakes, it is not clear that the Euler equations (or a simpler subset) 
should govern the far field behavior especially in the wake region. Because 
of this and the desire to provide a simple approach to this complicated 
problem, attention is focused in this work on modifying the physics of the 
problem. The modification is such that, if one uses time marching methods to 
obtain a steady state solution, the net result of the approach is to provide a 
different path of convergence to the steady state. 

Using index notation, the conservation of mass, momentum and energy 
equations can be written as 

^ + — U I = 0 (1) 
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where p is the density, Ui is the velocity component in the direction of x^, E 
and H are the total energy and enthalpy per unit mass, p is the pressure, ajj 
is the stress tensor, is the heat flux vector and Sjj is the Kronecker 
delta. Expressions for ajj and are needed to close the above system. In 
the traditional Navier-Stokes equations, the stress tensor and the heat flux 
vector have the fora 
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with p, K and K being the coefficients of shear viscosity, heat conductivity 
and bulk viscosity. The first order system of equations (1) - (5) cannot be 
hyperbolic because of the absence of time derivatives in equations (4) and 
(5). A necessary (but not sufficient) modification that may make the system 
hyperbolic is the addition of time dependent terms to equations (4) and (5), 
i.e. , change them to 
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where a and 0 are unknown parameters. It is seen from equations (7) and (8) 
that when the steady state is reached one will recover equations (4) and (5). 

The a and fi have to be chosen in such a way that the resulting system is 
hyperbolic. Rather than make aribtrary choices, attention was focused on the 
results of kinetic theory. As is well known, the stress tensor and the heat 
flux vector can be determined from the solution of the Baltzmann equation. 9 
The traditional Hilbert or Chapman-Enskog methods yield equations similar to 
equations (4) and (5). On the other hand, the thirteen moment method of 
Grad 10 yields equations similar to equations (7) and (8). By comparing 
equations (7) and (8) with the results of the thirteen moment method 
(equations (7.5-7) and (7.5-8) of Ref. 9), it is seen that a and f should be 
chosen as 


« = 1/P » # = l/c p p 

where c p is the specific heat at constant pressure. 


(9) 
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Using the above choices, the governing system of equations, (l)-(3), (7)- 


(8) can be written as 


»U 


A-+B- 


+ C 


dt dx ay 


au 

az 


= F(U) 


( 10 ) 


where 


U 





pB, 


o 


i y 



(id 


with A, B and C being appropriate matrices. In the most general case, some of 
the eigenvalues of A, B and C cannot be determined in cloaed form. It was 
necessary to use numerical computations to show that all of the eigenvalues 
were real. Thus, the resulting system is hyperbolic. 

The presence of source terms in the resulting formation renders the 
scheme som ew hat stiff. To overcome this problem a semi-implicit treatment, 
similar to that of Ref. 11, was employed. A finite volume, cell-centered, 
formulation and a four-step Runge-Kutta time stepping scheste, similar to that 
of Ref. 12, is used to obtain the solution. Local time stepping and residual 
smoothing 13 were employed to accelerate convergence to the steady state. 


Boundary Conditions 

In order to demonstrate the validity of this approach, it is necessary to 
compare its predictions with some of the well-known exact solutions. As an 
illustration, the flow past a semi-infinite plate (Blasius problem) is 
considered. Because the normal stresses and the axial heat flux are very 
small, a "thin layer" approximation in which e xx , Oyy and q x were assumed 
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negligible was developed. For this case, it was possible to derive analytical 
expressions for the eigenvalues of the matrices A and B. 

It may be recalled that the eigenvalues of A for the two-dimensional 
Euler equations are u, u, u*a with v replacing u for the eigenvalues of B. In 
the present case, where the thin layer approximation is used, the eigenvalues 
of A are 


u, u * a, 
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where a is the speed of sound. The matrix A can be written as 


A = P A P 


( 12 ) 


where A is a diagonal matrix whose elements are the eigenvalues of A, P is a 
matrix whose columns are the eigenvectors of A and P -1 is the inverse of P. 
The characteristic variables are given by 


W = P *u 
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where S is the entropy. 

At boundaries noraal to the x-axis, equation (10) aay be approximated by 


— + A — = F(U) (15) 

dt dx 


Multiplying equation (15) by P~ l , using eq. (13), and assuming that P ~ 1 is 
aprpoximately constant, one obtains 
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If a steady state solution is desired, the boundary conditions cure given by 


A — = P -1 F (17) 

dx 


Upwind differencing is used to express the derivatives at the boundary. 
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Results and Discussion 


Results for a flow past a flat plate at a Mach number of .3 and a 
Reynolds number of 5,000 are presented. Figure 1 shows a comparison of the 
dimensionless velocity in the axial direction, as a function of the Blasius 
coordinate with the Blasius solution. Similarly, Figure 2 compares the 
calculated Cf, the skin friction coefficient, as a function of Re x with the 
result obtained from the Blasius solution. Good agreement is indicated in 
both cases. Figure 3 is taken from Ref. 14 and is intended to compare our 
results with those of calculations based on the traditional Navier-Stokes 
equations. Again, good agreement is indicated. 

As a further validation of this approach we hope to include in the final 
paper the calculation of subsonic separated flow past a cylinder. 1 s The wake 
of the cylinder is dominated by a Karman-Vortex street. Because of lack of 
adequate boundary conditions, the problem proved to be difficult to simulate. 
The present approach, with boundary conditions provided by equation (16), and 
similar relations, is highly suited for the solution of this problem. 

In conclusion, we have developed a rather straightforward and general 
approach for determining inflow and outflow boundary conditions for the 
Navier-Stokes equations. It is hoped that this approach will make it possible 
to tackle more difficlut problems involving transition, separation and 
turbulance. 
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Figure 1. Comparison of the dimensionless velocity 
as a function of the Blasius coordinate. 
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Figure 2. Comparison of skin-firction distributions 
for laminar flat plate flow based on the 
Hyperbolic Navier-Stokes equations. 
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Figure 3. Comparison of skin-friction distributions 
for laminar flat plate flow based on the 
traditional Navier-Stokes equations. 


